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Abstract 

We present a consecutive derivation of mapping equations of motion for the one-dimensional 
classical hydrogenic atom in a monochromatic field of any frequency. We analyze this map in the 
case of high and low relative frequency of the field and transition from regular to chaotic behavior. 
We show that the map at aphelion is suitable for investigation of transition to chaotic behavior also 
in the low frequency field and even for adiabatic ionization when the strength of the external field 
is comparable with the Coulomb field. Moreover, the approximate analytical criterion (taking into 
account the electron's energy increase by the influence of the field) yields a threshold field strength 
quite close to the numerical results. We reveal that transition from adiabatic to chaotic ionization 
takes place when the ratio of the field frequency to the electron Kepler frequency approximately equals 
0.1. For the dynamics and ionization in a very low frequency field the Kepler map can be converted to 
a differential equation and solved analytically. The threshold field of the adiabatic ionization obtained 
from the map is only 1.5% lower than the exact field strength of static field ionization. 
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1. INTRODUCTION 

It is already the third decade when the highly excited hydrogen atom in a microwave 
field remains one of the simplest and most proper real system for experimental and theoretical 
investigation of classical and quantum chaos in the nonlinear systems strongly driven by external 
driving fields (see reviews [1-5] and references herein). For theoretical analysis of transition to 
stochastic behavior and ionization processes of atoms in microwave fields approximate mapping 
equations of motion, rather than differential equations, are most convenient. 

Such a two-dimensional map (for the scaled energy of the one-dimensional atom in a 
monochromatic field and for the relative phase of the field), later called Kepler map [3,13], 
has been obtained in Refs. [6,7] after an integration of equations of motion for one period of 
the intrinsic motion of the electron between two subsequent passings of the aphelion, the largest 
distance from the nucleus. This map greatly facilitates numerical investigation of dynamics and 
ionization process and allows even an analytical estimation of the threshold field strengths for 
the onset of chaos, the diffusion coefficient of the electron in energy space and other char- 
acteristics of the system [3-10]. Moreover, this map is closely related to the expressions of 
quasiclassical dipole matrix elements for high atomic states [11,12]. 

The Kepler map for relatively high frequencies of the field (recovered in [3,13] and sometimes 
represented for the number of absorbed photons) is relatively simple and is widely used for 
analysis of classical dynamics, as well as after the quantization - the quantum Kepler map for 
the quantum dynamics [3-5,9,10,13]. In the region of relatively low frequency of the field this 
map is sufficiently more complex, the threshold field strength for transition to chaotic behavior 
and ionization is considerably higher than that for transition to chaotic behavior in the medium 
and high frequency fields. 

Since derivation of the mapping equations of motion is based on the classical perturbation 
theory and the electron's energy change due to interaction with an external field during the 
period of motion in the Coulomb field depends on the initial condition, i.e., on the integration 
interval [6,7,11], complementary analysis of the applicability of the 'standard' Kepler map is 
necessary. 

It should be noted, that in derivation of the Kepler map one integrates not over the period 
of the external electromagnetic field but over the period of the electron intrinsic motion in 
the Coulomb field [6,7]. This results in some contradictions and difficulties. First, the period 
of integration and the obtained map depend on the energy of the electron which complicates 
the quantization problem of the Kepler map [3-7,13-17]. Second, the energy change of the 
electron during the period of the classical intrinsic motion due to interaction with the microwave 
field depends on the starting conditions, i.e., on the integration interval [7,11]. This makes it 
possibility to obtain another map, the Kepler map at perihelion, derived by the integration 
between two subsequent passings of the perihelion [7] or even a three-dimensional map for the 
two halves of the intrinsic period [11]. These maps stronger reveal the resonance structure of 
chaotic dynamics at low frequencies [6,7,11,15]. 

Moreover, Nauenberg [18] has presented a so-called canonical Kepler map which agrees 
with the results of Refs. [6,7,11] if taken at perihelion but not with the widely used Kepler 
map at aphelion. The expressions for the variables of this canonical Kepler map are, however, 
sufficiently complicated, not explicit and, therefore, they are not comfortable for analytical 
and even numerical analysis. So, for this reason corroboration of the standard Kepler map's 
applicability for the large range of parameters of the problem is significant as well. Note 
additionally, that in derivation and analysis of the maps in Ref. [7] some inaccuracy and 
misprints have appeared. 
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All these aspects indicate the need of additional analysis of the mapping equations of motion 
for a highly excited classical hydrogenic atom in a monochromatic field. In addition, transition 
from adiabatic to chaotic ionization mechanism in the a frequency field is of great interest (see, 
e.g. [19] and references herein). 

In this paper we present a consistent derivation of the mapping equations of motion for a 
one-dimensional classical atom in a monochromatic field of any frequency and analyze transition 
from regular to chaotic behavior and ionization process. From the fulfilled analysis we conclude 
that the map at aphelion and an approximate analytical criterion of the onset of chaos are 
suitable also in the low frequency region, even for adiabatic ionization, where the strength of 
the external field is comparable with the Coulomb field. Moreover, in this case the map can be 
transformed to a differential equation and solved analytically. 

2. MAPPING EQUATIONS OF MOTION 

The direct way of coupling the electromagnetic field to the electron Hamiltonian is through 
the A • P interaction, where A is the vector potential of the field and P is the generalized 
momentum of the electron. The Hamiltonian of the hydrogen atom in a linearly polarized field 
Fcos(ujt + ■&), with F, uj and $ being the field strength amplitude, field frequency and phase, 
respectively, in atomic units is 

H = ^ (p - ^ sin(cjt + tf^J -i. (1) 

The electron energy change due to interaction with the external field follows from the Hamil- 
tonian equations of motion [20] 

E = -r .Fcos(wf + #). (2) 

Note, that Eq. (2) is exact if r is obtained from equations of motion including influence of 
the electromagnetic field. Using parametric equations of motion in the Coulomb field we can 
calculate the change of the electron's energy in the classical perturbation theory approximation 
[6-8,11]. 

Measuring the time of the field action in the field periods one can introduce the scale 
transformation [16,18] where the scaled field strength and the scaled energy are F s = F/u 4 ^ 3 
and E s = E/u 2 ^ 3 , respectively. However, it is convenient [6-8,17] to introduce the positive 
scaled energy e = —2E S and the relative field strength F = Fn 4 , = F s /eq, with n being 
the initial effective principle quantum number, uq = (—2E ) 1 . The threshold values of the 
relative field strength F for the ionization onset depends weaker upon the initial effective 
principle quantum number n and the relative frequency of the field so = than the scaled 
field strength F s . 

We restrict our subsequent consideration to the one-dimensional model, which corresponds 
to states very extended along the electric field direction. Such classical one-dimensional model 
was first considered in Refs. [21] for the description of surface-state electrons, while a justifica- 
tion of the use of one-dimensional-like states for periodically driven hydrogen atoms appeared in 
[22]. Since that the one-dimensional model is widely used in theoretical analysis [2-11,13-17]. 

For the derivation of a map describing the motion of an electron in the superposition of the 
Coulomb and microwave fields we should integrate dynamical equations over some characteristic 
period of the system. The peculiarity of the system under consideration is that we are able to 
obtain explicit expressions for the change of the electron energy only for halves of the period 
and for the complete period, T = 2tt/(-2E) 3 / 2 = 2n/ uie ' 2 of the intrinsic electron motion in 
the Coulomb field [6-8,11] but not for the period of the external field. 
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Integration of Eq. (2) for motion between two subsequent passages at the aphelion (where 
x — and there is no energy exchange between the field and the atom) results in the change of 
the electron's energy [6,7,11] 

AE = — (irF/E) J' s (s) sintf. 

Here s = e~ z l 2 = u>/(— 2E) 3 ^ 2 = uj/Q is the relative frequency of the field, i.e., the ratio of the 
field frequency uj to the Kepler orbital frequency Q = (— 2E) 3 ^ 2 , and J' s (z) is the derivative of 
the Anger function with respect to the argument z. The derivative of the Anger function 



1 7r 

l' s (s) = — / sin [s (x — sinx)] sinrrda; 

7T J 



is a very simple analytical function which can be approximated quite well by some combination 
[7] of expansion in powers of s 



l + (5/24)s 2 . 



and of the asymptotic form 



-wf / \ b a smns 

where 

o i/3 2 2 / 3 

a = o/Q — — ~ 0.4473, b = — -—- ~ 0.41085. 

3 2 / 3 r (2/3) ' 3V3r (1/3) 

Introducing the scaled energy e = —2E/tu 2 / 3 and the relative field strength F = F/AEq we 
have 

As = -nF elh (e)sm$ (3) 

where e = —2E /uj 2 ^ 3 and 

h(e) = - £ J' s (s). (4) 
The change of the field phase d after the electron motion period in the Coulomb field is 

M = 2-kujT = 2vr/e 3/2 . (5) 

Defining the scaled energy and the phase before, £j,$j, and after, £j + i,i?j + i, passages of the 
electron of one intrinsic motion period we can introduce [23,24] a generating function G (£j+i, flj) 
of the map determined as 

£j = dG/dtfj, = dG/de j+1 . (6) 

In agreement with Eqs. (3) and (5) the generating function is (see also [3] for analogy) 

G (e j+1 ,i3j)= e j+1 $j - 47re j+ 1 1 /2 - 7vF e 2 h (e j+1 ) cos (7) 

and according to Eqs. (6) it generates the map 

f £ i+i = e i - *Foe 2 h (e j+1 ) sin^-, 

\ = -&j + 2Txje% - -KF Q e 2 r, (s j+1 ) cos^, 1 ' 
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Here 

^) = — 0) 

Note that the map (8) can be derived also without introduction of the generating function 
[6-8,11] but using the requirement of the area-preserving of the map (8) defined as 

* ( *r n° = L < 10 > 

It should also be noted that the map (14)-(19) in Ref. [7] is with the positive signs of terms 
in the right-hand site of Eq. (8) containing the field amplitudes F , i.e., it is derived for the 
reverse orientation of the atom with respect to the field orientation. Also note that a function 
sin^fc was inadvertently omitted from the right-hand side of Eq. (15) in [7]. 

The map (8) is the general mapping form of the classical equations of motion for the one- 
dimensional hydrogen atom in a microwave field derived in the classical perturbation theory 
approximation. Some analytical and numerical analysis of this map has already been done in 
Refs. [3,6-8]. Here we analyze different special cases of the map (8). 

3. HIGH FREQUENCY LIMIT 

For relatively high frequencies of the field, s ^> 1 (s > 2), theoretical analysis of the classical 
dynamics of the one-dimensional hydrogen atom in a microwave field is relatively simple. That 
is why the energy changes of the electron, (Ej +i — Ej) and {sj+\ — Ej), do not depend on the 
initial energy Ej and relative frequency s >> 1. Indeed, using the asymptotic form of the 
derivative of the Anger function, J' s (s) = b/s 2 ^ 3 , we have h(sj + i) = 46 = const., n (£j+i) = 
and, consequently, the following map 



Ej+i = Sj — 47r6F £Q sini?j, 
i+ i = $ J +2n/E 3 / + \ 



3/2 (11) 



Note, that scaled classical dynamics according to maps (8) and (11) depends only on single 
combination of the field parameters, i.e., on the scaled field strength F s = F el = F/u 4 ^ 3 (see 
also [16,17]). 

By the standard [23,24] linearization procedure, Ej = e + Aej, in the vicinity of the integer 
relative frequency (resonance), s = Eq 3 ^ 2 = m with m integer, the map (11) can be transformed 
to the standard (Chirikov) map 

Ij+i = Ij + Ksiw&j, 

>>j • i >>.i ■ 1j ■ i • (12) 

Here = -SttAej /e 5 /2 and K = 12n 2 bF /^. 

From the condition of the onset of classical chaos for the standard map, K > K c ~ 0.9816 
[1,23-25], we can, therefore, estimate the threshold field strength for chaotization of dynamics 
and ionization of the atom in the high frequency field 

F C = KJ (l2ir 2 bsl /3 ) ~ 0.02s 1/3 . (13) 

Sometimes one writes the map (11) for a variable iV = — l/2n 2 u>, which change gives the 
number of absorbed photons [3,13], 

N j+1 = Nj + 2vr (f/cj 5/3 ) sin^-, 
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■&j + 2nu (-2uN j+1 ) ' 



(14) 



We see that for such variables the dynamics of the system depends on two parameters: 
on the ratio F q = F/u 5 ^ 3 (in Refs. [16,17] F q = F/uj 5 ^ 3 was called the quantum scaled field 
strength) and on the field frequency oo. Map (14) is, therefore, not the most convenient one for 
analysis of the classical dynamics. 

In general there are, however, no essential difficulties in the theoretical analysis of classical 
nonlinear dynamics of the highly excited hydrogen atom in the microwave field of relative 
frequency So — — 0-5 when the field strength is lower or comparable with the threshold field 
strength for the onset of classical chaos, i.e., if the microwave field is considerably weaker than 
the characteristic Coulomb field. In such a case, the energy change of the electron during the 
period of intrinsic motion is relatively small and application of the classical perturbation theory 
for derivation of the Kepler map (8) is sufficiently correct. Further analysis of transition to 
chaotic behavior and of the ionization process can be based on the map (8) and for sq ~ 0.3-7-1.5 
results in the "impressive agreement" [5] between measured ionization curves and those obtained 
from the map (8) [5-10]. Even analytical estimation of the threshold field strengths based on 
this map is rather proper [6-8]. 

Considerably more complicated is the analysis of transition to stochastic motion and of 
ionization process in the region of low relative frequencies, sq < 0.3, [6-8,11]. 

4. LOW FREQUENCY LIMIT 

For the low relative frequencies of the microwave field, s<l, the map (8) can be simplified 
as well. Using expansion of the function 3' 3 (s) in powers of s, J' s (s) ~ s/2, for s <C 1 we have 
according to Eqs. (4) and (9) 



This map is slightly more complicated than map (11) for high frequencies, however, it can 
easily be analyzed numerically as well as analytically. Note first of all, that the energy change 
of the electron during the period of intrinsic motion (after one step of iteration), \£j+i — £j \ , 
is considerably smaller than the binding energy of the electron Ej ~ Eq if the field strength is 
lower or comparable with the threshold field strength, i.e., 2ttF (so/s^+i) — 2irF e 1 ^ 2 <C £o, 
or 27rF so <C 1 if F < ^o* — 0.13 and So <C 1. This indicates that the map (16) is probably 
suitable for description of dynamics even in the low frequency region where the field is relatively 
strong. 

In Figs 1 and 2 the results of the numerical analysis of maps (8) and (16) in the low frequency, 
s < 1, presented. We see that the threshold ionization field calculated from the maps 

approaches the static field ionization threshold Fq* ~ 0.13 when s ~~ * 0. This supports the 
presumption that the map (8) is valid even in the low frequency limit where the strength of the 
driving field is of the order of the Coulomb field. 

4.1. Adiabatic ionization 

For low frequencies, 2ns = 2n/s 3 ^ 2 <C 1, according to the second equation of map (16) the 
change of the angle i? after one step of iteration is small. As it was noticed above, the energy 



h (e j+1 ) = 
V (£j+i) = 

Consequently map (8) transforms to the form 




(15) 




(16) 
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change is also relatively small. Therefore, we can transform the difference equations (16) to 
differential equations of the form 



de 27T£gF . 

d7 = "^7^ sm ^' 

^ = % + ^r cos ^ (17) 

Dividing second equation of the system (17) by the first one we obtain one differential equation 

d(cos^) e 5cos-# 

The analytical solution of Eq. (18) with the initial condition e = eo when ■& = $ is 

cos^ = z 5 cos^o-22 4 (l-2)/F , z = Je/e~ . (19) 

Eq. (16) describes the motion of the system in e and d variables, i.e., represents the 
functional interdependence between two dynamical variables. 

Let us analyze Eqs. (18) and (19) in more detail. For relatively low values of F , i.e., for 

F < |z 4 = | (j^j , the right-hand side of Eq. (18) is positive for all phases •&. Therefore, 
cos-# and e decrease and increase simultaneously and, according to Eq. (16), there is a motion 
in the whole interval [0, 2n] of the angle i?. For Fq > |z 4 , however, the increase of the angle •& 
in the interval -j- tt goes to decrease of d at d ~ ir. This results in fast decrease of e and to 
the ionization process (see also Fig. 1). It is easy to understand from analysis of Eq. (19) that 
the minimal value of Fq for such a motion (resulting in ionization) corresponds to i?o = and 
$ = 7r. This value of F is very close to the maximal value of F resulting from motion in the 
whole interval [0, 2ir] of i.e., the maximum of the expression 

F = 2z (1 - z) I (l + z b ) . (20) 

This maximum is at z = zo, where zo is a solution of the equation z 5 + 5z — 4 = 0, being 
zq — 0.75193. The critical value of the relative field strength, therefore, is Fq = 2zq/5 = 0.1279 
which is only 1.5% lower than the adiabatic ionization threshold Fq 1 = 2 10 / (37r) 4 = 0.1298 
[7,21]. According to our numerical analysis, if so < 0.05 the electron remains bounded and the 
dynamics is regular for F < 0.13 while ionization takes place for F > 0.131 (see also Fig. 1 and 
2). These results are very close as well to analytical conclusions. Note, that some decrease of 
the threshold field strength values Fq with decreasing of s was observed for s < 0.1. Dynamics 
and classical ionization at such frequencies are, however, essentially adiabatic. 
4.2. Chaotic ionization 

For higher relative frequencies, s > 0.1, ionization process is due to chaotic dynamics of 
the highly excited electron of the hydrogenic atom in a microwave field. There are different 
criterions for estimation of the parameters when the dynamics of the nonlinear system becomes 
chaotic. For analysis of transition to chaotic behavior of the motion described by maps (8), 
(11) and (16) the most proper, to our mind, is the criterion related with the randomization of 
the phases [24] 

5# j+1 



K = max 



1 



> I. (21) 
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Here max means the maximum with respect to the phase -&j and variation of the phase with 
respect to the phase -&j means the full variation including dependence of on dj through 
the variable Ej+i in Eqs. (8), (11) and (16). 

Applying criterion (21) to the general map (8) we obtain the threshold field strength 



£ 7/2 



If £ ~ e Eq. (22) yields the result 



F ° C -12vr%^( S )- (22) 



F c =(l27r 2 S J' s ( S )) _1 (23) 



which for s ^> 1 coincides with Eq. (13). 

For more precise evaluation of the critical field strengths we should take into account the 
change (increase) of the electron's energy due to the influence of the electromagnetic field. For 
higher relative frequency s or lower scaled energy Ej the threshold ionization field is lower. 
Therefore, if the scaled energy Sj decreases as a result of relatively regular dynamics in a not 
very strong microwave field, then the lower field strength is sufficient for transition to chaotic 
dynamics. For high frequencies such change of the energy is relatively small. Nevertheless it 
reveals some resonance structure in the field-atom interaction. In the low frequency limit the 
energy change is more essential. Now consider it in detail. 

As it was shown above, maximal decrease of the scaled energy Ej is for the angle $j ~ n and 
it can be evaluated from Eq. (20). Taking this into account we have from Eq. (16) according 
to criterion (21) the expression for the threshold relative field strength 



e 5 zl° 



F ° _ 6vr%2 " 67r 2 s 2 (24) 

where z c is the solution of Eq. (20) with F = F C . Eq. (20) can be solved approximately 
expanding in powers of Ffi. The result of such an expansion is 

zl° ~ 1 - 10F C + 30 (F C ) 2 - 73.6 (F C ) 3 (25) 

where the last term in the right-hand side of Eq. (25) is from the requirement of the exact 
maximal value z c = Zq = 0.75193 for the static threshold field strength Fq = 0.1279. 

For evaluation of the threshold field for transition to chaotic behavior in the low frequency 
field we should thus solve the system of equations (24) and (25). For 0.09 < s < 0.5 expressions 
(24) and (25) give an ionization threshold field very close to the numerical results (see Fig. 2). 
For frequencies lower than sq ~ 0.09 ionization is adiabatic, because for so low frequencies the 
adiabatic ionization threshold field, Fq = 2z^/5 = 0.1279, is lower than the phase randomization 
field evaluated according to Eqs. (24) and (25). The adiabatic ionization, therefore, occurs in 
such a case earlier than the chaotization of the dynamics. Note, that numerical results reveal 
transition from adiabatic to chaotic ionization at relative frequency s — 0.1 (scaled energy 
E ~ 4.3) as well. 

At higher frequencies ionization is due to chaotic dynamics while transition to chaotic be- 
havior can be evaluated from the approximate criterion (21) taking into account the electron's 
energy change by the influence of the electromagnetic field. For frequencies higher than s — 0.5 
we should use a more exact expression than J' s (s) ~ s/2 for the derivative of the Anger function, 
i.e., Eq. (22). 
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5. CONCLUDING REMARKS 

From the analysis given in this study we can conclude that the map at aphelion (8) is 
suitable for investigation of regular and stochastic classical dynamics, transition to chaotic 
behavior and ionization of Rydberg atoms in high, medium and low frequency fields, even for 
adiabatic ionization when the strength of the external field is comparable with the averaged 
Coulomb field. For such a purpose it is unnecessary to use the map at perihelion, the map 
for two halves of the intrinsic period or the canonical Kepler map [7,11,18]. Moreover, the 
approximate criterion (21) for transition to chaotic behavior yields a threshold field strength 
very close to the numerical results if we take into account increase of the electron energy by 
influence of the electromagnetic field. Transition from adiabatic to chaotic ionization of the 
classical hydrogenic atom in a monochromatic field takes place at a relative field frequency 
s ~ 0.1. 

Furthermore, the Kepler map and some generalizations of it (for two- and multi-frequency 
[9,10,26] or some other fields, e.g., circular polarized microwave field, for three-dimensional 
atoms and other modifications) are and may be more widely used for analysis of different effects 
of classical and quantum chaos in driven nonlinear systems [19,27,28]. Note also the attempts 
to derive and use similar maps in astronomy for analysis of chaotic dynamics of comets and 
other astronomical bodies [29]. It turns out, however, that in such a case, generalization of 
the Kepler map for nonharmonic perturbations and for motion in three-dimensional space, is 
necessary. 
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Caption for the figures 



Fig. 1. Trajectories of the map (8) for different initial conditions, e , 9 , and different relative 
field strength Fq. The pictures in the left-hand side correspond to the regular quasiperiodic 
motion while those in the right-hand side represent ionization process for a little stronger field. 
At e ~ 4.3, i.e., s ~ 0.11 a transition from the adiabatic to the chaotic ionization mechanism 
takes place. 

Fig. 2. Relative threshold field strength for the onset of ionization from numerical analysis 
of the maps (8) and (16) and according to the approximate criterion (24) - (25). 
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Fig. 1. Trajectories of the map (8) for different initial conditions, £ , 0o, and different relative 
field strength Fo- The pictures in the left-hand side correspond to the regular quasiperiodic 
motion while those in the right-hand side represent ionization process for a little stronger field. 
At £o=4.3, i.e., a so=0.11 a transition from the adiabatic to the chaotic ionization mechanis 
takes place. 



